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We study the discovery reach for Non-Standard Interactions (NSIs) in a neutrino factory exper- 
iment. After giving a theoretical, but model- independent, overview of the most relevant classes of 
NSIs, we present detailed numerical results for some of them. Our simulations take into account 
matter effects, uncertainties in the neutrino oscillation parameters, systematical errors, parameter 
correlations, and degeneracies. We perform scans of the parameter space, and show that a neutrino 
factory has excellent prospects of detecting NSIs originating from new physics at around 1 TeV, 
which is a scale favored by many extensions of the standard model. It will also turn out that the 
discovery reach depends strongly on the standard and non-standard CP violating phases in the 
Lagrangian. 
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1. INTRODUCTION 

Huge efforts are currently undertaken to design new 
long-baseline neutrino experiments to precisely measure 
the three-flavor oscillation parameters, in particular the 
yet unknown mixing angle 0i3, the CP violating phase 
(5cPj and the sign of the atmospheric mass squared dif- 
ference Ato|]^. However, the excellent accuracy with 
which the planned setups can measure the oscillation 
probabilities, will also allow for the detection of new 
sub-leading effects, such as mixing with sterile neutri- 
nos, a non-unitary Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix, neutrino decay, the decoherence ef- 
fect, CPT violation, or mass-varying neutrinos. Further- 
more, many extensions of the standard model predict 
new, effective four-Fermi interactions involving neutri- 
nos, on which we will focus in this paper. General phe- 
nomenological studies of these non-standard interactions 
(NSIs) have been conducted in |2|, |3j. A, ^ , ,61, 0, Q, 
and specific models are discussed in [3, After a 

work by Grossman which pointed out the impor- 
tance of NSIs for neutrino oscillation experiments, many 
authors have investigated their impact in the context 
of solar neutrinos fl2. [isl [Tst . atmospheric neutri- 
nos [l^, [13, m, Il9l . 12(1 l2 l|. conventional and upgraded 
neutrino beamsj2|j2lj24^,'2l,[2^^ neutrino facto- 
ries d, m. Hi, [mlull, Mi , beta beams [H , super- 
nova neutrinos l36l . l37j| cosmological relic neutrinos [3^ , 
e+e" colliders |39|, neutrino-electron scattering [43|, and 
neutrino-nucleus scattering [4ll. l4^ . Existing experimen- 
tal bounds are presented in |43j . 

In this article, we will discuss in particular the dis- 
covery potential of a neutrino factory [H, [H, 
m, [43, [H^l , which is currently the most advanced tech- 
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nology discussed in neutrino physics, and would have a 
precision of 0(1-0.1)% on some of the oscillation proba- 
bilities. We will show that this translates into a sensitiv- 
ity to NSIs originating from scales of up to several TeV. 
We will focus only on non-standard interactions (NSIs) 
which violate lepton flavor. Existence of such interactions 
would typically induce not only effects in the neutrino 
sector, but also charged lepton flavor violating processes 
like fj, 3e. However, in the charged lepton sector, the 
signal is proportional to the square of the non-standard 
coupling, while in an oscillation experiment, interference 
between the standard and non-standard amplitudes will 
also induce terms which are linear in the coupling con- 
stant and can therefore be expected to be easier to detect. 



A long baseline neutrino oscillation experiment con- 
sists of three stages: beam production (source), beam 
propagation through the Earth, and neutrino detection 
at the far site. Here, we are going to consider NSIs which 
modify only one of these aspects at a time. In principle, 
there may also exist combined effects of several new pro- 
cesses, but these will be suppressed by higher powers of 
the small coupling constants. 



The paper is organized as follows: In Sec.[21 we will de- 
scribe the NSIs analytically, but in a model-independent 
way. Afterwards, in Sec. [31 we will present detailed nu- 
merical results on non-standard modifications to the neu- 
trino production and propagation amplitudes. We have 
performed sophisticated simulations with a modified ver- 
sion of the GLoBES software [Sll, , taking into account 
systematical errors and correlations between all standard 
and non-standard oscillation parameters. Some of these 
correlations will turn out to be very strong, so our final 
results on the NSI discovery reach of a neutrino factory 
will strongly depend on the true parameter values. We 
will summarize our results in Sec. S] and draw some con- 
clusions. 
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2. NON-STANDARD INTERACTIONS IN 
NEUTRINO OSCILLATIONS 

In the context of neutrino factory experiments, one 
usually considers the "golden oscillation channel" v^. — > 
Vu, , the "silver channel" i^e ^ i^t, the "platinum channel" 



Ve, and the disappearance channel i/^j 



(see 



e.g. [47[ and references therein). Of these, the golden 
channel is most important for the discovery of small ef- 
fects such as 013-oscillations, CP violation, but also non- 
standard interactions, because it is an appearance chan- 
nel, and because it is technically more easily accessible 
than the silver and platinum channels. Therefore, we 
will focus on the golden channel and consider only NSIs 
influencing the corresponding process chain 



Osc. 



(1) 



NSIs can modify the production, oscillation, and detec- 
tion of neutrinos, so that the following alternative pro- 
cesses to Eq. ^ can occur: 



NSI 



No osc. 



No osc. 
NSI 



NSI 



(2) 
(3) 

(4) 



These processes are illustrated diagrammatically in 

Fig.m 

Since the initial and final states are the same in Eq. ([1]) 
and Eqs. ^ to ffl jinterference can occur on the level of 
the amplitudes 0, S, i, i, S, HH, |3l|. This will enhance 
the magnitude of the new effects compared to scenarios 
where the NSIs are added non-coherently. 

Non-standard interactions are typically generated by 
new physics at very high energy scales, so for a neutrino 
factory operating in the low-energy regime, they can be 
expressed as model-independent four-Fermi interactions: 

+ E ^^r/^K7^(l-7'K}{/7p(lT7')/} 

f—e,u,d 



Gf 
V2 



At 



(5) 



Here, Gp denotes the Fermi constant, Ve and are the 

neutrino fields, e and n are the charged lepton fields, 

and M, d are quark fields. Finally, the magnitude of the 

NSIs is parameterized by e^^ for effects in the neutrino 

source, by e™'-''^ for non-standard matter effects on the 

oscillation, and by e'if for modifications to the detection 

[ — h 

process. In contrast to a previous study [31[, we allow 
the e parameters to be complex. Then, Eq. ^ implies 
that the e parameters for anti-neutrino processes are the 
complex conjugates of those for neutrino processes. 



The expected magnitude of the e parameters can be 
estimated as follows t29|] : If we assume the non-standard 
interactions to arise at a scale Mnsi, the effective vertices 
in Eq. ([5]) will be suppressed by in the same 

way as the standard weak interactions are suppressed by 
Therefore we expect 



^NSl' 



(6) 



The NSIs in the beam source, given by the first line 
of Eq. ^ , effectively promote the initial neutrino state 
from a pure flavor eigenstate into the mixed state 



(7) 



Similarly, the detector will project out the mixed state 



EuEd 



(8) 



rather than the pure flavor eigenstate {I'^.l- The kinemat- 
ics of a neutrino factory experiment is such that in most 
regions of the phase space, -Ee,/j,ti,d ^ Tne,^,u,d- This re- 
mains true even if we take into account that Eg and as 
well as Eu and Ed are not independent, and if we allow E^ 
and Ed to be slightly off-shell due to QCD effects in the 
nucleus. We do not consider these hadronic effects here, 
but we remark that they will generally be soft compared 
to the primary neutrino interaction. Since interference 
of standard and non-standard processes can only occur 
if all initial and final state particles have the same chi- 
rahty, it follows that the {V - A){V + A) type NSIs are 
suppressed by the hehcity factors O {m^rrif,/ E^E^) resp. 
O {rriumd/ EuEd). Therefore, we will take 



(9) 
(10) 



in the following. 

Note that the e parameters do not necessarily form 
unitary matrices, so that the source and detection states, 
in general, do not form complete sets of basis vectors in 
the Hilbert space: 



E 



l^i'^X^i'^l^i- (11) 



However, we do require the PMNS matrix to be unitary, 
so that the standard mass and flavor eigenstates at least 
form a basis of the subspace of states participating in 
oscillations. Thus, the neutrino propagation does not vi- 
olate unitarity, while the production and detection pro- 
cesses may do so. With these assumptions, neutrino os- 
cillations can be described as usual by a hermitian 3x3 
Hamiltonian 13111. that contains, however, an extra term 
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H^si due to the second line of Eq. . Thus we can write 



H ^ H; 



so 



/o 



so 



2E 



(C/t) 



2E 



2E 



^NSi(e"^) 



2E ' 



(12) 

(13) 
(14) 



Here, i?so contains the standard osciUations (SO), and 
a is the effective matter potential, which we assume to 
be constant in the following. The effective NSI coupling 
e™ is related to the parameters e^^^'-^ from Eq. ^ by 
the formula 



m,e- 



+ 361^-). (15) 



This relation can be understood if we assume the num- 
bers of protons, neutrons, and electrons to be the same in 
the Earth matter and consider only the effect of valence 
quarks. Furthermore, we have made use of the fact that 
the spin and momentum average of the Earth is zero, so 
that only the components of the vector interactions in 
the second line of Eq. ([5]) are relevant, and contribute 
equally. 

In principle, -ffNSi can also contain other non-zero en- 
tries besides e™ . In combination with standard oscilla- 
tions, these can lead to process chains like 



No osc. 

i^e > l^T 

NSI 



Osc. 



(16) 



In the remainder of this section, we will neglect such con- 
tributions for conciseness, but we will exemplarily con- 
sider effects proportional to e™ in our numerical analysis 
in Sec. 13.21 A systematic study of non-standard Hamil- 
tonians like Eq. is given in 3 • 

The amplitude of the flavor transition Va ^ i^p can be 
calculated from the propagation Hamiltonian Eq. 
byi 

(17) 

Here, the standard oscillation amplitude 6*30 reads 

5so = e-'^^°^, (18) 

and the amplitude induced by the non-standard matter 
effects, S'nsIi is given by the perturbative expansion 



{'S'NSi(e™ )}/3q - ^ [Sso] 



^ ) ■ya 



(19) 



In our case, the (e,^) and (/^,e) elements of 
{-f^NSi(e"^)}/3Q are non-zero, so the golden-channel fla- 
vor transition i^e f/i can occur even in the absence of 
standard oscillations. 



^ Here we regard the NSI parameters CejT*' to be small pertur- 
bations. In |26|I the authors pointed out that, from the current 
experimental limits, the NSIs might even dominate over the os- 
cillation effect in a — ► Ur oscillation experiment. In such a 
situation, this perturbative expansion would no longer be valid. 
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If, as a last step, we replace the initial and final states 
in Eq. by the modified states from Eqs. Q and (H]), 
we obtain the transition probability up to first order in 
the e parameters as 



(20) 



|(5'so)/.eP 

+ 2|e™|Re[(5so);e{^NSi(e™)U 

(21) 



2|e:jRe 



2k?^|Re 



0(e^ 



The zeroth order term represents the standard oscilla- 
tion probability, while the first order terms contain the 
contributions from the diff'erent types of NSIs. 

In this work, we are interested in the discovery poten- 
tial for non-standard effects, i.e. in the prospects of iden- 
tifying the tiny NSI contribution on the large standard 
oscillation background. If only the terms proportional to 
are present, this can be achieved by exploiting the 
different spectral structure of the signal and background 
events [23j. If we expand the oscillation amplitudes up 
to first order in l/E, we find that {Sso)fie ^ ^/E, while 

- 1. 



according to Eq. {5nsi(C)} 



('S'so)a 



Hence, the first (standard) term in Eq. (PT|) behaves as 
while the second (non-standard) term is propor- 
tional to l/E. The situation is quite different for effects 
proportional to e™, since for these, the non-standard 
terms {S'NSi(e^)};ie will contain a factor (5'so)mt ~ 1/^: 
so their lowest order energy dependence is identical to 
that of the standard oscillations. Therefore, the discov- 
ery reach of a neutrino factory for e™ will be worse than 
that for e™ . 

For non-standard effects parameterized by e^^ and ef^, 
we can read off from Eq. (|2ip that we are again in a 
favorable situation, since {Sso)*^e{Sso) tifi ~ ^/E and 



3. DETECTING NON-STANDARD 
INTERACTIONS IN A NEUTRINO FACTORY 

To obtain reliable estimates for the prospects of dis- 
covering non-standard interactions in a neutrino factory, 
we have performed detailed numerical simulations with 
a modified version of the GLoBES software [Bll, H^]. We 
use a neutrino factory setup based on NuFact2 from [ssj . 
with a parent muon energy of 50 GeV and a baseline of 
3000 km. The total running time is 8 years (4 years in 
the neutrino mode, 4 years in the anti-neutrino mode), 
and the number of stored muons per year is 1.066 • 10^^. 
The detector is a 50 kt magnetized iron calorimeter, and 
the cross sections are based on [HI, . Both the wrong- 
sign muon appearance channel ( "golden channel" ) and 
the muon disappearance channel are taken into account. 





i/e appearance 


ff, disappearance 


Signal 


2.5% 


20% 


Background 


2.5% 


20% 



Table I: Systematical flux normalization uncertainties in our 
neutrino factory setup NuFact2. 



We have incorporated the backgrounds due to neutral 
current events and muon charge misidentification. 

We quantify the performance of an experiment by in- 
troducing the discovery reach for non-standard interac- 
tions, which is defined as the minimal magnitude of the 
e parameters, for which the expected experimental data 
is no longer consistent with a standard oscillation fit. 

Following the statistical procedure described in the ap- 
pendix of [5J|, we define the following function ^ 



Y = mm 
A 



where Nj, 



channel bin |^^^ ^^true^ ^true) _ j^^^ (A, e ^ 0)|^ 



7Vy(A'™%e'™<=) 

+ Priors, (22) 



denotes the number of events in the z-th 
energy bin for the oscillation channel j, the vector 
A = (6*12, 6*13, 6*23, (Sep, ^^21, Am§2, a, fe) contains the 
standard oscillation parameters, the Mikheyev-Smirnov- 
Wolfenstein (MSW) potential a, and the systematical bi- 
ases b, and e represents the non-standard parameters. 
The index j runs over the — > v^, and t'^ — *■ v^, channels 
and over the corresponding anti-neutrino processes. For 
the "true" parameters used to calculate the simulated 
data, we adopt the following numerical values [56| : 



sin^ 261*™° 


= 0.83, 


sin^ 261*™° 


= 1.0, 


sin^ 261*™° 


= 0.01, 


(Am^i)'™'= 


= 8.2 X 10"^ eV^ 


(Amii)*™° 


= 2.5 X 10"^ eV^ 



In the fit, we marginalize over all standard oscillation 
parameters and over the systematical biases, but keep 
the non-standard parameters fixed at 0. The prior terms 
implement external input from other experiments and 



have the form {x 



')^/(T^, where x stands for any 



oscillation parameter or systematical bias, and ax is the 
corresponding externally given uncertainty. We assume 
012 and Am^x to be known to within 10% from solar and 



^ In the actual implementation, we assume the events to follow 
the Poisson distribution. However, for illustrative purposes, it is 
sufficient to consider the more compact approximative Gaussian 
expression. 
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reactor experiments [56*], and include a standard devia- 
tion of 5% for the MSW potential a. All other oscilla- 
tions parameters are assumed to be unconstrained since 
the neutrino factory itself has an excellent sensitivity to 
them. The systematical uncertainties are summarized in 
Tab.m 

For compactness, our discussion will focus on non- 
standard interactions induced by e™ , e™ , and e^^, but of 
course, one could also derive similar results for the other 
possible terms. In particular, NSIs in the detector can be 
expected to have similar effects to those in the source. 

Furthermore, we will always assume a normal mass hi- 
erarchy, both for the simulated data and for the fit. The 
main influence of the inverted hierarchy is to shift the at- 
mospheric MSW resonance to the anti-neutrino channel, 
which is in general less important for the overall sensitiv- 
ity of the experiment because of the smaller anti-neutrino 
cross section. However, one can easily see that the dis- 
covery reach for non-standard interactions is robust with 
respect to the presence or absence of the MSW resonance: 
The main effect of the resonance is to enhance |S'so| in 
Eq. (HI]). Therefore, if it is effective, the signal term, 
which is proportional to jS'gg'S'NSili is enhanced. At the 
same time, however, also the standard oscillation back- 
ground proportional to |S'soP will become larger. These 
two opposing effects cancel each other, as can be seen 
from the expression (P^ : If we assume A*''"'' = A, the 
background terms drop out in the numerator, but not 
in the denominator. Since, however, the numerator con- 
tains an extra square, we obtain ~ I'S'gg'S'NSiP/l'S'soP, 
i.e. the standard oscillation contributions cancel, and the 
expression is unaffected by their MSW enhancement. We 
have verified numerically, that our results would hardly 
be affected by using the inverted hierarchy for the data 
and the fit, even if we included parameter correlations 
and higher order terms. 



3.1. Effects proportional to 

We will first concentrate on non-standard effects pro- 
portional to e™ , and assume all other e parameters to 
vanish. Fig. O shows the NSI discovery reach at 3cr as 
a function of the true values of 6cp and arg[e™], and 
for three different values of |(e^)'™°|- Since, according 
to Eq. (P^ . arg[e™] and |e™ | are fixed at zero in the 
fit, while all other parameters are marginalized over, the 
contours are based on the assumption of two degrees of 
freedom.'^ If the parameters lie in the white regions of 
the plots, non-standard interactions can be established 



^ One might argue that the leading term in the oscillation prob- 
ability depends on the parameter combination arg[€™] -|- 5qp 
(see Eq. II23I I below), so it may be justified to use only 1 d.o.f. 
However, since sub-leading contributions are not completely neg- 
ligible, we take 2 d.o.f. to be conservative. 



at the 3cr level, while in the dark areas, the sensitivity is 
less than la. It is obvious that for larger |(ej|!^ |, the 
white regions of good sensitivity become larger. 

The characteristic band structure in Fig. [2] reveals that 
there are strong correlations between 6cp and arg[e™]. 
To understand these correlations analytically, we note 
that the leading NSI signal term is proportional to 
cos(arg[e™] + 6cp) (Hj. Therefore, the contribution of 
the i^fi channel to the function from Eq. (PH) 

becomes approximately 

2 {0\r ■ l(ee"^)'™1 • cos (arg[(6™ + 5*,™°))' 

(23) 

and is thus expected to be roughly constant along the 
lines of constant arg[(e™ )*''"°] -t-J^p". This behavior can 
be nicely seen in Fig. [2l 

Note that correlations do not only limit the discovery 
reach for non-standard interactions, but can also compli- 
cate the measurement of the standard oscillation param- 
eters j32j . 

Comparing the three plots in Fig. [21 we find that for 
|(e™)'''"''| ~ 6 • 10"'', the first white islands appear, 
i.e. there are some parameter combinations for which 
the non-standard effects can be discovered at Sc. For 
|(-gm ^truoj > 4 . ^^le x^ valucs are above 3cr in the 

whole parameter space, i.e. a 3(t discovery is always pos- 
sible, independent of arg[(e™ )*''"°] and Sq^'^. According 
to Eq. ([S]), these numbers translate into a sensitivity to 
mass scales of up to A/nsi ~ 1 — 3 TeV. 

These values reappear as the top and bottom edges of 
the foremost (green) bars in the middle part of Fig. [31 
The bars stacked below them (light blue) show how the 
discovery reach would improve if all standard oscillation 
parameters and the MSW potential were known with in- 
finite precision, and the hindmost (blue) bars have been 
calculated have been calculated under the additional as- 
sumption that systematical errors are not present. The 
plot shows that the discovery reach depends crucially on 
the true values of arg[e™] and Jcp, while systematical 
errors and the correlations with the fit parameters have 
only moderate impact. 

Comparing the results for different values of 6113, we 
find that the achievable sensitivity for the most favorable 
combination of phase parameters (top edges) remains 
roughly unchanged as ^13 decreases, while that for the 
most problematic parameters (bottom edges) becomes 
slightly better. On the one hand, smaller ^13 means 
smaller (S'so)Ate, so all terms in Eq. (|2ip will decrease. 
On the other hand, the standard oscillation background, 
which is given by |(S'so)/ieP and therefore proportional 
to 0^3, will decrease faster than the non-standard term, 
which is linear in Oi^ [2^. This makes it easier to disen- 
tangle signal and background, and especially when corre- 
lations between the standard and non-standard param- 
eters are taken into account, this improved background 
suppression seems to overcompensate the smaller signal. 
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|(e")tt"'=|=6-10-'* 
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Figure 2: Contour plots of the function defined in Eq. ^ in the 5j^^J'=-arg[(e™ plane for |(e™ )*''"°| = 6 x IQ-* (top), 
1 X 10""^ (center), and 4 x 10""^ (bottom). The value of sin^ 2^13 was taken to be 10~^ in all cases. The standard oscillation 
parameters and the matter potential were marginalized over. 

Finally, let us compare the sensitivities predicted by parabolic ^, this translates into a 3cr bound < 1.4 • 10~^. 
our simulations with existing bounds from charged lepton 

flavor violating processes. In [i^, Davidson et al. con- 

strain |e™ | to be smaller than ~ 8 • 10~* at the 90% con- 
fidence level. Assuming the respective function to be We are aware, that this extrapolation is problematic, since real- 
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(e^^)'™^ discovery reach (3cr) for NuFact-II 

Full marginalization GLoBES 

SO parameters fixed 

Systematics off 




Figure 3; Limitations to the discovery reach for |e™ | arising from systematical errors and from parameter correlations. The top 
edges of the bars indicate the values of for which there exists some combination of arg[(e™ )*'^"°] and ^cp"^, which yields 

a sensitivity better than 3a; the bottom edges show how large the NSIs need to be in order to be detected at this confidence 
level for all possible values of arg[(e™ )''^"'^] and ^cp"^. The foremost (green) bars were obtained with the full analysis procedure 
discussed in the text, while for the intermediate (light blue) bars, the marginalization over standard oscillation parameters 
(including the MSW potential) was omitted, and for the hindmost (blue) bars, systematical errors were also switched off. 



Comparing this number with Fig. [31 we find that part of 
the parameter space accessible by the neutrino factory 
is already ruled out, but, depending on the phase cor- 
relation, our setup may still have a significant discovery 
potential. This is particularly interesting if we note that 
the present 90% C.L. bound, |e™ | < 8- 10"'*, comes actu- 
ally from three independent bounds on the coherent for- 
ward scattering on up-quarks (< 8 • 10^'^), down-quarks 
(< 8 • 10"''), and electrons (< 5 • 10"''). The former 
two have been derived from jj, e conversion in nuclei, 
while the latter stems from the constraints on — s- 3e. 
To improve the overall bound on |e™ | significantly, all 
three components would have to be improved. Since no 
experiment searching for /i — > 3e is being designed at the 
moment, a neutrino factory seems to be the most realistic 
future option for studying |e™ | in a model-independent 
way. Of course, if the LHC should find evidence for one 
specific class of models, the present bound might become 
much stronger already within the next few years. 



istic functions can be far from parabolic. However, it should 
give a useful order of magnitude estimate. 



3.2. Effects proportional to 

Let us now turn to non-standard effects proportional to 
e™ , which are introduced in analogy to Eqs. - (fH)) . 
It can be read off from Fig. SI that the sensitivity of a neu- 
trino factory to these effects is almost two orders of mag- 
nitude worse than that to : Only for |e™ | > 3 • 10'^ 
discovery can be guaranteed. This can be understood 
from our discussion in Sec. [H which shows that the en- 
ergy dependence of standard and non-standard effects is 
the same, so the effect of e™ can easily be absorbed into 
A. This also explains why fixing the standard oscilla- 
tion parameters improves the sensitivity by one order of 
magnitude. 

As in the case of |e™ |, the ^13 dependence in Fig. [3] is 
weak. Actually, it could be expected to be larger here, 
because for the large values of |e™ | required for discovery, 
not only the interference term between standard oscilla- 
tions and NSI will contribute to the oscillation probabil- 
ity, but also the pure NSI term proportional to |e™p. 
For large ^13, both terms are comparable in magnitude, 
while for small ^13, but still large |e™ |, the pure NSI term 
dominates. Therefore, in the latter case, the qualitative 
behavior of the transition probability could be expected 
to be simpler and thus more easily absorbed into Sep- 
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(e^^)'*^"^ discovery reach (3cr) for NuFact-II 

Full marginalization GLoBES 

SO parameters fixed 

Systematics off 




sin^20l'^'= = lO-^ sin^20^i'^" = lO-^ sin^26'l'^" = 10-^ 



Figure 4: Limitations to the discovery reach for |e™ | arising from systematical errors and from parameter correlations. The 
color-coding is the same as in Fig. O 



However, we see from Fig. U that this effect is not very 
pronounced. 

The present bounds on |e™ | are of [HI, [Hi, so the 
neutrino factory could break new ground, independent of 
the phase correlations. Note that in [3l|, the authors pre- 
dict an even better performance for the neutrino factory. 
However, they employ a completely different experimen- 
tal setup with a baseline of only 732 km, and use a differ- 
ent analysis technique. Note also that it has been pointed 
out in [ij, [s^l that a measurement of the Weinberg angle 
by neutrino scattering in the near detector of a neutrino 
factory could improve the limit on | e™ | independently to 
well below 0.1. 

3.3. Effects proportional to e^/j 

If the non-standard interactions do not affect neutrino 
oscillations, but rather the production process, we expect 
from Sec.[2]that the sensitivity will again be excellent be- 
cause the standard and non-standard terms have differ- 
ent energy dependence. Indeed, Fig. [S] shows that effects 
with jcg^l ^ 10~^ might be detected, and that detection 
can be guaranteed for |eg^| > 10"^. 

This discovery reach is at least one order of magnitude 
better than the model-independent bound of 0(10"^) 
coming from universality considerations in lepton de- 
cays [29j . Let us, however, remark, that present model- 
dependent bounds on e^^ can be already stronger than 



the sensitivity of the neutrino factory. 



4. CONCLUSIONS 

We have investigated the prospects of a search for non- 
standard neutrino interactions in a neutrino factory ex- 
periment. We have discussed several different contri- 
butions that can arise in the effective Lagrangian, and 
have pointed out that these can be distinguished from 
the standard oscillations by their characteristic energy 
dependence. We have performed careful numerical sim- 
ulations of a neutrino factory experiment to estimate its 
discovery potential for e™ , e™, and e^^. It turned out 
that there is a strong correlation between arg[e™] and 
6cp , so that the discovery reach for e™ depends strongly 
on the true values of these parameters: For certain com- 
binations, a discovery of the non-standard interactions is 
possible for |e™ | < 10~^, while for less favorable scenar- 
ios, |e™ I ~ 10~^ is required. Since the present bounds on 
|e™ I are already of O(10~'^), the discovery potential will 
crucially depend on the phase correlations. Vice-versa, 
however, a combination of neutrino factory results with 
limits from other experiments might provide additional 
constraints on the phases. The sensitivity to |e™ | is more 
than one order of magnitude worse than that to |e™ | due 
to the less favorable energy dependence of this effect. 
However, since present bounds on |e™ | are very weak, a 
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Figure 5: Limitations to the discovery reach for |e^^| arising from systematical errors and from parameter correlations. The 
color-coding is the same as in Fig. O 



neutrino factory could achieve a significant improvement 
here. Finally, the sensitivity to \el^\ ranges between 10"'^ 
and 10~^, which is at at least one order of magnitude 
better than present model-independent bounds. Thus, 
our simulations show that a neutrino factory is an excel- 
lent tool for detecting new physics in the neutrino sector. 
However, reversing the argument, this also means that 
possible non-standard interactions have to be taken into 
account when analyzing the data of such an experiment. 



Acknowledgments 



We would like to thank P. Huber and W. Winter for 
useful discussions. This work was in part supported by 
the Transregio Sonderforschungsbereich TR27 "Neutri- 
nos and Beyond" der Deutschen Forschungsgemeinschaft. 



[1] L. Wolfenstein, Phys. Rev. D17, 2369 (1978). [11] 

[2] J. W. F. Valle, Phys. Lett. B199, 432 (1987). 

[3] M. M. Guzzo, A. Masiero, and S. T. Petcov, Phys. Lett. [12] 

B260, 154 (1991). 
[4] E. Roulet, Phys. Rev. D44, 935 (1991). 

[5] S. Bergmann, Y. Grossman, and E. Nardi, Phys. Rev. [13] 

D60, 093008 (1999), hep-ph/9903517. 
[6] T. Hattori, T. Hasuike, and S. Wakaizumi, Prog. Theor. [14] 

Phys. 114, 439 (2005), hep-ph/0210138. 
[7] M. Garbutt and B. H. J. McKellar (2003), hep- [15] 

ph/0308111. 

[8] M. Blennow, T. Ohlsson, and W. Winter (2005), hep- [16] 
ph/0508175. 

[9] A. De Gouvea, G. F. Giudice, A. Strumia, and K. Tobe, [17] 
Nucl. Phys. B623, 395 (2002), hep-ph/0107156. 
[10] T. Ota and J. Sato, Phys. Rev. D71, 096004 (2005), hep- [18] 
ph/0502124. 



Y. Grossman, Phys. Lett. B359, 141 (1995), hep- 
ph/9507344. 

S. Bergmann, M. M. Guzzo, P. C. de Holanda, P. L 
Krastev, and H. Nunokawa, Phys. Rev. D62, 073001 
(2000), hep-ph/0004049. 

Z. Berezhiani, R. S. Raghavan, and A. Rossi, Nucl. Phys. 
B638, 62 (2002), hep-ph/0111138. 

A. Friedland, C. Lunardini, and C. Pena-Garay, Phys. 
Lett. B594, 347 (2004), hep-ph/0402266. 
O. G. Miranda, M. A. Tortola, and J. W. F. Valle, JHEP 
10, 008 (2006), hep-ph/0406280. 

M. C. Gonzalez-Garcia et al, Phys. Rev. Lett. 82, 3202 
(1999), hep-ph/9809531. 

S. Bergmann, Y. Grossman, and D. M. Pierce, Phys. Rev. 
D61, 053005 (2000), hep-ph/9909390. 
N. Fornengo, M. Maltoni, R. T. Bayo, and J. W. F. Valle, 
Phys. Rev. D65, 013010 (2002), hep-ph/0108043. 



10 



[19] M. C. Gonzalcz-Garcia and M. Maltoni, Phys. Rev. D70, 

033010 (2004), hcp-ph/0404085. 
[20] A. Friodland, C. Lunardini, and M. Maltoni, Phys. Rev. 

D70, 111301 (2004), hep-pli/0408264. 
[21] A. Friedland and C. Lunardini, Phys. Rev. D72, 053009 

(2005) , hep-ph/0506143. 

[22] S. Bergmann and Y. Grossman, Phys. Rev. D59, 093005 

(1999), hcp-ph/9809524. 
[23] T. Ota, J. Sato, and N.-a. Yamashita, Phys. Rev. D65, 

093015 (2002), hcp-ph/0112329. 
[24] T. Ota and J. Sato, Phys. Lett. B545, 367 (2002), hep- 

ph/0202145. 

[25] M. Honda, N. Okamura, and T. Takeuchi (2006), hep- 
ph/0603268. 

[26] N. Kitazawa, H. Sugiyama, and O. Yasuda (2006), hep- 
ph/0606013. 

[27] A. Friedland and C. Lunardini, Phys. Rev. D74, 033012 

(2006) , hep-ph/0606101. 

[28] M. Blennow, T. Ohlsson, and J. Skrotzki (2007), hep- 
ph/0702059. 

[29] M. C. Gonzalcz-Garcia, Y. Grossman, A. Gusso, and 

Y. Nil, Phys. Rev. D64, 096006 (2001), hep-ph/0105159. 
[30] P. Huber and J. W. F. Valle, Phys. Lett. B523, 151 

(2001), hep-ph/0108193. 
[31] A. M. Gago, M. M. Guzzo, H. Nunokawa, W. J. C. Teves, 

and R. Zukanovich Funchal, Phys. Rev. D64, 073003 

(2001), hcp-ph/0105196. 
[32] P. Huber, T. Schwctz, and ,J. W. F. Valle, Phys. Rev. 

D66, 013006 (2002), hcp-ph/0202048. 
[33] M. Campanelli and A. Romanino, Phys. Rev. D66, 

113001 (2002), hep-ph/0207350. 
[34] A. Bueno, M. Campanelli, M. Laveder, J. Rico, and 

A. Rubbia, JHEP 06, 032 (2001), hcp-ph/0010308. 
[35] R. Adhikari, S. K. Agarwalla, and A. Raychaudhuri, 

Phys. Lett. B642, 111 (2006), hep-ph/0608034. 
[36] G. L. FogU, E. Lisi, A. Mirizzi, and D. Montanino, Phys. 

Rev. D66, 013009 (2002), hep-ph/0202269. 
[37] H. Duan, G. M. Puller, J. Carlson, and Y.-Z. Qian (2006), 

astro-ph/0608050. 



[38] G. Mangano et al., Nucl. Phys. B756, 100 (2006), hep- 

ph/0607267. 

[39] Z. Berezhiani and A. Rossi, Phys. Lett. B535, 207 
(2002), hep-ph/0111137. 

[40] J. Barranco, O. G. Miranda, C. A. Moura, and J. W. F. 
Valle, Phys. Rev. D73, 113001 (2006), hep-ph/0512195. 

[41] J. Barranco, O. G. Miranda, and T. L Rashba, JHEP 12, 
021 (2005), hep-ph/0508299. 

[42] J. Barranco, O. G. Miranda, and T. I. Rashba (2007), 
hep-ph/0702175. 

[43] S. Davidson, C. Pena-Gaxay, N. Rius, and A. Santamaria, 
JHEP 03, Oil (2003), hep-ph/0302093. 

[44] S. Geer, Phys. Rev. D57, 6989 (1998), hep-ph/9712290. 

[45] C. Albright et al. (2000), hep-ex/0008064. 

[46] M. ApoUonio et al. (2002), hcp-ph/0210192. 

[47] P. Huber, M. Lindner, M. Rolinec, and W. Winter, Phys. 
Rev. D74, 073003 (2006), hep-ph/0606119. 

[48] P. Huber, T. Schwctz, and J. W. F. Valle, Phys. Rev. 
Lett. 88, 101804 (2002), hep-ph/0111224. 

[49] M. Freund, P. Huber, and M. Lindner, Nucl. Phys. B615, 
331 (2001), hep-ph/0105071. 

[50] M. Freund, P. Huber, and M. Lindner, Nucl. Phys. B585, 
105 (2000), hep-ph/0004085. 

[51] P. Huber, M. Lindner, and W. Winter, 
Comput. Phys. Commun. 167, 195 (2005), 
http : //www . ph . turn . de/~ globes , hep-ph/0407333. 

[52] P. Huber, J. Kopp, M. Lindner, M. Rolinec, and W. Win- 
ter (2007), hep-ph/0701187. 

[53] P. Huber, M. Lindner, and W. Winter, Nucl. Phys. 
B645, 3 (2002), hep-ph/0204352. 

[54] M. D. Messier, Evidence for neutrino mass from observa- 
tions of atmospheric neutrinos with Super-Kamiokande, 
Ph.D. thesis, Boston University (1999), uML99-23965. 

[55] E. A. Paschos and J. Y. Yu, Phys. Rev. D65, 033002 
(2002), hep-ph/0107261. 

[56] M. Maltoni, T. Schwetz, M. A. Tortola, and J. W. F. 
Valle, New J. Phys. 6, 122 (2004), hep-ph/0405172. 

[57] M. L. Mangano et al. (2001), hep-ph/0105155. 



